Critical temperature for quenching of pair correlations 
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The level density at low spin in the ^^^■^^■^Dy and ^^^-^'^^Yb nuclei has been extracted hom primary 
7 rays. The nuclear heat capacity is deduced within the framework of the canonical ensemble. The 
heat capacity exhibits an S-formed shape as a function of temperature, which is interpreted as a 
fingerprint of the phase transition from a strongly correlated to an uncorrelated phase. The critical 
temperature for the quenching of pair correlations is found at = 0.50(4) MeV. 
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The thermodynamical properties of nuclei deviate from 
infinite systems. While the quenching of pairing in su- 
perconductors is well described as a function of temper- 
ature, the nucleus represents a finite many body system 
characterized by large fluctuations in the thermodynamic 
observables. A long-standing problem in experimental 
nuclear physics has been to observe the transition from 
strongly paired states, at around T = 0, to unpaired 
states at higher temperatures. 

In nuclear theory, the pairing gap parameter A can 
be studied as function of temperature using the BCS 
gap equations 0J^. From this simple model the gap de- 
creases monotonically to zero at a critical temperature of 
Tc ~ 0.5 A. However, if particle number is projected out 
|^,D, the decrease is significantly delayed. The predicted 
decrease of pair correlations takes place over several MeV 
of excitation energy ||] . Recently ||] , we reported struc- 
tures in the level densities in the 1-7 MeV region, which 
are probably due to the breaking of nucleon pairs and a 
gradual decrease of pair correlations. 

Experimental data on the quenching of pair correla- 
tions are important as a test for nuclear theories. Within 
finite temperature BCS and RPA models, level density 
and specific heat are calculated for e.g. ^*Ni within 
the shell model Monte Carlo method (SMMC) |^,|| one 
is now able to estimate level densities H in heavy nuclei 
1^ up to high excitation energies. 

The subject of this Letter is to report on the observa- 
tion of the gradual transition from strongly paired states 
to unpaired states in rare earth nuclei at low spin. The 
canonical heat capacity is used as a thermometer. Since 
only particles at the Fermi surface contribute to this 
quantity, it is very sensitive to phase transitions. It has 
been demonstrated from SMMC calculations in the Fe 
region ]ri|-p^, that breaking of only one nucleon pair 
increases the heat capacity significantly. 

Recently |lj,|l^ , we presented a method for extracting 
level density and 7 strength function from measured 7 
ray spectra. Since the 7 decay half lives are long, typi- 
cally 10~^^-10~"'^^ s, the method should essentially give 
observables from a thermalized system p^ . The spin 
window is typically 2-8 h and the excitation energy res- 
olution is 0.3 MeV. 



The experiments were carried out with 45 MeV ^He 
projectiles from the MC-35 cyclotron at the University 
of Oslo. The experimental data were recorded with the 
CACTUS multidetector array ||l^ using the (^IIe,a7) re- 
action on I62,i63j)y g^j-^^ i72a73Yi3 self-supporting targets. 
The beam time was two weeks for each target. The 
charged ejectiles were detected with eight particle tele- 
scopes placed at an angle of 45° relative to the beam 
direction. Each telescope comprises one Si Ai? front and 
one Si(Li) E end detector with thicknesses of 140 and 
3000 /im, respectively. An array of 28 5"x5" Nal(Tl) 7 
detectors with a total efficiency of ^15% surrounded the 
target and particle detectors. 

From the reaction kinematics the measured a particle 
energy can be transformed to excitation energy E. Thus, 
each coincident 7 ray can be assigned a 7 cascade orig- 
inating from a specific excitation energy. The data are 
sorted into a matrix of {E,E^) energy pairs. At each 
excitation energy E the Nal 7 ray spectra are unfolded 
JlSf , and this matrix is used to extract the primary 7 ray 
matrix, with the well established subtraction technique 
ofRef. 

The resulting matrix P{E,Ej), which describes the 
primary 7 spectra obtained at an initial excitation energy 
E, is factorized according to the Brink-Axel hypothesis 
by P{E, E^) (X p{E-E-y)a{E^). The level density 
p(E)a.nd the 7 strength function (j{E^) are determined 
by a least fit to P. Since the fit yields an infinitely 
large number of equally good solutions, which can be ob- 
tained by transforming one arbitrary solution by 
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we have to determine the parameters A, B and a by 
comparing the p and a functions to known data. The A 
and a parameters are fitted to reproduce the number of 
known levels in the vicinity of the ground state and 
the neutron resonance spacing at the neutron binding 
energy 



Figure Dj shows the extracted level densities and 7 
strength functions for the ^^^'^^"^Dy and I'^i'i'^^yi, nuclei. 
The data for the even nuclei are published recently 
and are included to be compared to odd nuclei. Also the 



1 



7 strength functions are given for all four nuclei. Since 
other experimental data on the 7 strength function in 
the energy interval 1.5-6 MeV is sparse, we could not 
fix the parameter B\ and therefore we give the functions 
in arbitrary units. However, when fitting the functions 
with the expression Ci?" we obtain n ~ 3.9, 4.5 for the 
Dy and Yb isotopes respectively, which one would expect 
from the tale of the giant dipole resonance |l4) . In the 
following, only the level density will be discussed. 
The partition function in the canonical ensemble 



Z{T) = Y^p{E^)e-^-l^ 



(3) 



n=0 



is determined by the measured level density of accessible 
states p{En) in the present nuclear reaction. Strictly, the 
sum should run from zero to infinity. In this work we 
calculate Z for temperatures up to T = 1 MeV. Assum- 
ing a Bethe-like level density expression |2^] , the average 
excitation energy in the canonical ensemble 



(EiT)) = J2 Enp{En)e~^-^^ 



(4) 



gives roughly (E) ~ aT^ with a standard deviation of 
as ^ T\/2aT where a is the level density parameter. 
Using Eq. (|) requires that the level density should be 
known up to {E) + Sct^, typically 40 MeV. However, the 
experimental level densities of Fig. |l| only cover the ex- 
citation region up close to the neutron binding energy of 
about 6 and 8 MeV for odd and even mass nuclei, re- 
spectively. For higher energies it is reasonable to assume 
Fermi gas properties, since single particles are excited 
into the continuum region with high level density. There- 
fore, due to lack of experimental data, the level density is 
extrapolated to higher energies by the shifted Fermi gas 
model expression p^] 



and the averaging made in Eq. (Q) gives a smooth temper- 
ature dependence of Cy(r). A corresponding Cv{{E)) 
may also be derived in the canonical ensemble. 

The deduced heat capacities for the I6i462j-)y 
i7i,i72y]-, nuclei are shown in Fig. ||. All four nuclei ex- 
hibit similarly S-shaped Cv(T)-curves with a local maxi- 
mum relative to the Fermi gas estimate at Tc « 0.5 MeV. 
The S-shaped curve is interpreted as a fingerprint of a 
phase transition in a finite system from a phase with 
strong pairing correlations to a phase without such cor- 
relations. Due to the strong smoothing introduced by the 
transformation to the canonical ensemble, we do not ex- 
pect to see discrete transitions between the various quasi- 
particle regimes, but only the transition where all pairing 
correlations are quenched as a whole. In the right pan- 
els of Fig. ^, we see that Cv{{E)) has an excess in the 
heat capacity distributed over a broad region of excita- 
tion energy and is not giving a clear signal for quenching 
of pairing correlations at a certain energy ||^. 

In order to extract a critical temperature for the 
quenching of pairing correlations from our data, we have 
to be careful, not to depend too much on the extrapola- 
tion of p. An inspection of Fig. |l| shows that the level 
density is roughly composed of two components as pro- 
posed by Gilbert and Cameron (i) a low energetic 
part; approximately a straight line in the log plot, and 
(ii) a high energetic part including the theoretical Fermi 
gas extrapolation; a slower growing function. For illus- 
tration, we construct a simple level density formula com- 
posed of a constant temperature level density part with 
T as temperature parameter, and a Fermi gas expression 



p{E)^ 



r/eicp{E/T) 
E-^/^ exp(2> 



for £' < e 
for E > e 



(7) 



where t] — e I exp(2Y^— e/r) accounts for continuity 
at the energy E = e. If we also require the slopes to be 
equal at e, the level density parameter a is restricted to 



12V(rT776ai/2f/3/2 ^1/3' 



(5) 



where U is the shifted energy and A is the mass num- 
ber. For shift and level density parameters a, we use 
von Egidy's parameterization of this expression |2§|| . The 
expression had to be normalized by a factor / in order to 
match the neutron resonance spacing data. The factors 
are 1.2, 0.94, 0.4 and 0.6 for I6i462j3y ^^^^ i7i472Yb re- 
spectively. The solid lines in Fig. |l| show how the expres- 
sion extrapolates our experimental level density curves. 

The extraction of the microcanonical heat capacity 
Cv (E) gives large fluctuations which are difficult to inter- 
pret g. Therefore, the heat capacity Cy(T) is calculated 
within the canonical ensemble, where T is a fixed input 
value in the theory, and a more appropriate parameter. 
The heat capacity is then given by 



Cv{T) 



d{E) 
dT ' 



2V^ 



(8) 



Figure || shows the heat capacity evaluated in the 
canonical ensemble with the level density function of 
Eq. (|) and = 1.7 MeV'^ The left hand part sim- 
ulates a pure Fermi gas description, i.e. the case e = 0, 
assuming a level density parameter a = 20 MeV~^. One 
can see that a pure Fermi gas does not give rise to the 
characteristic S-shape of the heat capacity as in Fig. ||. 
The right hand part simulates the experiments, where 
e = 5 MeV and a fulfills Eq. (|), i.e. again a = 20 MeY'^. 
The characteristic S-shape emerges. 

Therefore, our method to find Tc relies on the assump- 
tion that the lower energetic part of the level density can 
be approximately described by a constant temperature 
level density. Calculating {E{T)) and Cy(T) within the 
canonical ensemble for an exponential level density gives 



(6) T-^ = {E{T))-^ + ' 



and 
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Cv{T) = {\~T/t)-\ (9) 

Thus, plotting as function of {E(T))~^ one can de- 
termine T from Fig. ^. The quantity r is then identified 
with the critical temperature Tc, since Cy(T) accord- 
ing to Eq. (^) exhibits a pole at r and the analogy with 
the definition of Tx in the theory of superfluids becomes 
evident. The CviT) curve of Eq. (|) with = t us- 
ing the extracted critical temperatures for the four nu- 
clei is shown as dashed-dottcd lines in Fig. ||. This sim- 
ple analytical expression with only one parameter Tc fits 
the experimental data up to temperatures of ^--^0.4 MeV. 
The critical temperature itself is marked by the vertical 
lines. The extracted Tc is rather close to the estimate of 
Tc represented by the arrows. The extracted values are 
Tc = 0.52, 0.49, 0.50 and 0.49 MeV for the I6i,i62j3y 
i7i,i72y]3 nuclei respectively, which are somehow delayed 
compared to a degenerated BCS model with Tc — 0.5 A 
yielding Tc - 0.48, 0.46, 0.41 and 0.38 MeV for the re- 
spective nuclei, where A is calculated from neutron sep- 
aration energies p^ . 

We will now discuss how sensitive the extracted criti- 
cal temperature is with respect to the extrapolation, and 
we will give an estimate of the uncertainty of the ex- 
tracted critical temperatures. In the fit in Fig. ^ we use 
only energies from (E) ~ 0.5-2 MeV. This corresponds 
to energies in the level density curves up to i?„ ~6 MeV 
according to Eq. (jj). Also in Fig. ^, the interval where 
Eq. (|) fits the experimental data is T ~0-0.4 MeV. This 
corresponds to energies in the level density curves up to 
En ~8 MeV. Thus, the extracted critical temperature 
does only depend weakly on the actual extrapolation of 
p curve. Also the S-shape of the Cy(r) curve depends 
only on the fact, that the nuclear level density develops 
somehow as a Fermi gas expression at energies somewhere 
above the neutron binding energy. However, the actual 
values of the CviT) curve above T ~0.5 MeV do depend 
on the specific extrapolation chosen. 

With respect to Tc, the extrapolation is only impor- 
tant for determining the parameters A and especially a 
of Eq. (||) . We can therefore estimate the error of Tc by 

(10) 

where AE is the energy difference between the upper and 
lower energy where A and a are determined. Aa is the 
uncertainty of the level density parameter, AU is the en- 
ergy difference between the neutron binding energy and 
the upper point where A and a are determined, [/„ is 
the (shifted) neutron binding energy. D and AD are the 
neutron resonance spacing and its error. The errors of 
5% are added in order to account for the two fitting pro- 
cedures, one fitting A and a the other fitting Tc, both 
with an uncertainty of some 5%. Using a very conser- 
vative estimate of a « 17.5(6.0) MeV"\ C/„ < 8 MeV, 
AU < 3 MeV, D/AD < 0.2 and AT: > 4.5 MeV, we 



obtain ATc < 0.04 MeV. This yields a maximum error 
of Tc of some 8%. It is also important to notice, that 
due to the strong smoothing in Eq. (^ , the errors of the 
experimental level density curves are negligible in our 
calculation. 

In conclusion, we have seen a fingerprint of a phase 
transition in a finite system for the quenching of pair- 
ing correlations as a whole, given by the S-shape of the 
canonical heat capacity curves in rare earth nuclei. For 
the first time the critical temperature Tc at which pair 
correlations in rare earth nuclei are quenched, has been 
extracted from experimental data. The reminiscence of 
the quenching process is distributed over a 6 MeV broad 
excitation energy region, which is difficult to observe and 
interpret in the microcanonical ensemble. Simple argu- 
ments show that the peak in the heat capacity arises from 
two components in the level density; a constant temper- 
ature like part and a Fermi gas like part. It would be 
very interesting to compare our results with SMMC cal- 
culations performed for a narrow spin window. 
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FIG. 1. Experimental level density (points in left panels) and 7 strength function (right panels) for ^''^'^^^Dy (upper part) 
and "I'l^'^Yb (lower part). The error bars show the statistical uncertainties. The solid lines are extrapolations based on a 
shifted Fermi gas model (see text). The isolated point at the neutron binding energy is obtained from neutron resonance spacing 
data. 
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FIG. 2. Semi-experimental lieat capacity as function of temperature (left panels) and energy (E) (right panels) in the 
canonical ensemble for ^^^'^^^Dy and ^^^'"'^'^■^ Yb. The dashed lines describe the approximate Fermi gas heat capacity. The 
arrows indicate the first local maxima of the experimental curve relative to the Fermi gas estimates. The dashed dotted lines 
describe estimates according to Eq. (bl) where t is set equal to the critical temperature Tc. Tc is indicated by the vertical lines. 
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FIG. 4. The critical temperature is deduced by fitting a straight Une with slope 1 to the data points between the arrows. 
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